Charge order in Kondo lattice systems 
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Charge order is a commonly observed phenomenon in strongly correlated materials. However, 
most theories are based on a repulsive inter-site Coulomb interaction in order to explain charge order. 
We here show that only due to local interactions, charge order is favorable in heavy fermion systems 
at quarter filling. The driving force is the Kondo effect, which leads to a non-linear dependence on 
the filling of the lattice site for the energy gain. We use dynamical mean field theory combined with 
numerical renormalization group to demonstrate the existence of charge order at quarter filling for 
a Kondo lattice model. 

PACS numbers: 71.10.Fd;71.27.+a;75.25.Dk 



I. INTRODUCTION 

Heavy fcrmions are one of the most studied materials 
in condensed matter physics. This is due to a large num- 
ber of intriguing phenomena ranging from quantum crit- 
icality to unconventional superconductivity, which can 
be found in these systems Many of these phenom- 
ena are caused by a competition between the Kondo ef- 
fect and long range order, as visualized in the Doniach 
phase diagram^ Most of the long-range ordered phases, 
which have been studied theoretically so far in these ma- 
terials, are magnetically ordered phases like antifcrro- 
magnetism or ferromagnetism. However, besides these 
phases charge order has also been found experimentally 
in heavy fermion systems, e.g. in Rsli^Siio (R = Lu 
or Yb) .— These compounds exhibit a charge density wave 
(CDW) at high temperatures, and a cooperation between 
the CDW and a magnetic or superconducting state at 
much lower temperatures. Charge order can be naturally 
expected in systems at quarter filling for strong enough 
repulsive inter-site Coulomb interaction. This repulsive 
inter-site Coulomb interaction will lead to a situation, 
where it is less favorable for electrons to occupy neigh- 
boring sites equally than forming a CDW. The forma- 
tion of a CDW thus minimizes the effect of the repulsive 
inter-site interaction. A similar effect can be observed by 
introducing a repulsive c- /-electron interaction in heavy 
fermion models. 

We will show here that it is possible to stabilize a 
CDW in the Kondo lattice model without any inter-site 
Coulomb interaction. As for many other interesting ef- 
fects in heavy fermion systems, the driving force is the 
Kondo effect. Due to the exponential dependence of the 
energy, which is gained by the Kondo effect, it can be en- 
ergetically favorable for the system to form a CDW. The 
possibility of a CDW in the Kondo lattice model was 
already proposed by Hirsch— in a strong-coupling expan- 
sion, who showed that perturbation theory creates an ef- 
fective repulsive inter-site interaction, which may lead to 
a CDW phase. Recently, Otsuki et al.— have shown that 



the static susceptibility for a CDW diverges at quarter 
filling for small coupling strengths in the Kondo lattice 
model, indicating the existence of a CDW. Furthermore, 
phases with cooperation of spin and charge order have 
been reported for heavy fermion models on geometri- 
cally frustrated lattices, and are known as partial Kondo 
screened phases^— In these cases, some sites in the unit 
cell form a magnetically ordered state, while others stay 
magnetically disordered, but may have a different occu- 
pation number. 

The purpose of this paper is to show that a pure CDW 
can be obtained at quarter filling without frustration un- 
der the assumption of a bipartite lattice, and determine 
the physical properties of this phase. Furthermore, we 
will propose an explanation for the stability of a CDW 
in terms of the Kondo effect. 

For studying the CDW in heavy fermions, we perform 
calculations for the Kondo lattice model. The Kondo 
lattice model describes a lattice of local moments which 
are coupled to itinerant electrons. The Hamiltonian can 
be written as 4,13,14 
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where the first term describes electron hopping on a lat- 
tice and the second term the spin-spin interaction be- 
tween a localized moment S and the spin of a conduc- 
tion electron, s*j = cj p B pp iCi p i (with Pauli matrices a). 
Throughout the paper we will assume an antiferromag- 
netic coupling, J > 0, between the conduction electrons 
and the localized spins, which is appropriate to describe 
heavy fermion systems. 

For solving the Kondo lattice model we use the single 
site dynamical mean field theory (DMFT)J^^ DMFT 
relates the lattice model to a quantum impurity model, 
which must be solved self-consistently. As a consequence 
of this mapping onto an impurity model, the momentum 
dependence of the self-energy in the obtained solution 
is neglected. This approximation becomes exact in in- 
finite dimensions. However, comparison to experiments 



2 



has shown that for three dimensional systems DMFT rep- 
resents a good approximation, which is able to explain at 
least qualitatively many effects seen in real experiments. 
To solve the effective quantum impurity model, we use 
the numerical renormalization group (NRG) j 17,18 NRG 
is able to accurately calculate spectral functions for arbi- 
trary parameters and wide range of temperatures , 19 i 20 
and is also able to resolve exponentially small energy 
scales, which might emerge in quantum impurity mod- 
els. 

This paper is organized as follows: In the next section, 
we will show the existence of the CDW state at quar- 
ter filling by analyzing solutions without magnetic order. 
This will allow us to directly study the properties of the 
charge ordered state. In the third section, we allow mag- 
netic and charge order at quarter filling simultaneously, 
and analyze the cooperation/competition between these 
two forms of long range order. 
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Figure 1: (Color online) Charge order parameter for a param- 
agnetic state at quarter filling for different temperatures and 
coupling strengths J. Lines are a guide for the eye. 



II. PARAMAGNETIC 
CHARGE-DENSITY- WAVE AT QUARTER 
FILLING 

The phase diagram of the Kondo lattice model within 
DMFT has been investigated by a number of authors 
beforei^— Around quarter filling, most previous stud- 
ies showed the existence of a magnetically long range or- 
dered state. However, recent quantum Monte Carlo cal- 
culations of static susceptibilities by Otsuki et al.— have 
shown the divergence of the CDW susceptibility around 
quarter filling. 

In order to go beyond the study of static susceptibil- 
ities and to investigate the properties of the charge or- 
dered state at quarter filling, we perform calculations for 
an infinite dimensional Bethc lattice with bandwidth W 
and DOS p{w) = 8/(ttVF 2 ) y/W 2 /4 - uj 2 . Furthermore, 
we divide the lattice into two sublattices, from now on 
called A- and B-sublattice. In order to calculate the 
CDW state, we initialize these two sublattices with differ- 
ent chemical potentials in the first DMFT iteration, and 
then iteratively solve the DMFT self-consistency equa- 
tion, using the same chemical potential for both lattice 
sites. Furthermore, breaking the spin symmetry by ap- 
plying a magnetic field in the first DMFT iteration will 
allow to study combinations of charge order and mag- 
netism. Preserving the spin symmetry will lead to a 
paramagnetic state. 

In this section, we neglect magnetically ordered states. 
We show in Fig. [1] that a CDW exists for the Kondo 
lattice model even without inter-site repulsive Coulomb 
interaction. All solutions exhibiting charge order, which 
are found in this study, are located in the vicinity of 
quarter filling, (tia) + (ub) = 1- In Fig. [TJwe show the 
charge order parameter An = \{tia) — (it-b)\ for differ- 
ent temperatures and coupling strengths. Furthermore, 
a phase diagram plotting the amplitude of the CDW over 
temperature and coupling strength is shown in Fig. [5] 
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Figure 2: (Color online) Magnitude of the charge order pa- 
rameter for different temperatures and coupling strengths. 



Increasing the coupling strength J, the transition tem- 
perature as well as the charge order parameter itself in- 
creases. The magnitude of the charge order reaches a 
maximum of \ua ~ ns\ ~ 0.23 for J ps 0.5VF. Thus, 
the less occupied lattice sites are even at T — filled by 
(n) « 0.38. At weak coupling strengths, J/W < 0.3, the 
transition temperature increases rapidly with increasing 
coupling strength. Such a rapid change in the transition 
temperature could be explained by the exponential de- 
pendence of the Kondo effect on the coupling strength. 
Further evidence that the Kondo effect is essential for 
stabilizing charge order at quarter filling is that we find 
no sign of charge order for the ferromagnetically coupled 
Kondo lattice model, J < 0, for which the Kondo effect 
is not present. 

At a coupling strength J w (1 /2) W, the behavior 
changes dramatically. The transition temperature and 
the amplitude of the order parameter begin to decrease. 
More importantly, charge order seems to disappear at 
T = 0. For coupling strengths 0.5W < J < 0.65W we 
find a reentrant behavior (homogeneous - charge order - 
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homogeneous) for decreasing the temperature. Eventu- 
ally, for J ps 0.65VF charge order vanishes completely for 
all temperatures. 

To clarify the behavior at strong coupling strengths in 
more detail, we show in Fig. [3] the order parameter at 
fixed temperature in the lower panel and the region of 
chemical potential, for which a charge ordered state can 
be found at quarter filling, in the upper panel. It is evi- 
dent from Fig. [3] that the width of the parameter region, 
in which the CDW can be found, decreases with increas- 
ing coupling strength. The amplitude of the charge order 
within this region is only slightly changed when increas- 
ing the coupling strength. However, at a critical coupling 
strength the parameter region vanishes completely. This 
behavior indicates a first order phase transition at T = 0. 
The corresponding jump in the order parameter is shown 
in the lower panel of Fig. [3] At high temperatures, on the 
other hand, the order parameter decreases gradually to 
zero, indicating a second order phase transition. Thus, at 
the critical interaction strength and finite temperatures, 
the phase transition changes from first to second order. 

For J/W > 0.55, we do not find a CDW at T = 0. 
Only at temperatures T/W > 0.004, we can stabilize 
charge order. As we will show in Fig. [6l the charge or- 
dered state possesses a gap at the Fermi energy. The 
region of chemical potential, in which the charge ordered 
state can be stabilized at quarter filling, is supposed to 
be related to the gap width, at least at T = 0. This is be- 
cause the state does not change when changing the chem- 
ical potential within this region. Therefore, the shrink- 
ing of the parameter region corresponds to a decrease of 
the gap width near the critical coupling strength. Fur- 
thermore, we see that the gap is closed for J/W = 0.6 
for T/W < 0.004. Even if there would be a tiny re- 
gion of chemical potential stabilizing charge order for 
J/W ps 0.6 at T = 0, the gap width of the state would be 
tiny and this state would be unstable towards perturba- 
tions. Thus, we believe that a reentrant behavior occurs 
for strong coupling strength in the vicinity of a first order 
phase transition. 

Comparing to experiments, reasonable interaction 
strengths are J/W < 0.3. If we assume a bandwidth 
of W = leV in our calculations, the resulting Kondo 
temperature as well as the transition temperature of the 
charge order for these interaction strengths are T < 
WOK, which are realistic temperatures for /-electron sys- 
tems. 

J.E. HirscbJi already pointed out that there exists a 
nearest-neighbor repulsive density-density interaction in 
second-order strong coupling expansion, which for weak 
coupling might be sufficient to drive the system into a 
CDW. However, the applicability of a strong coupling 
theory to explain the existence of the CDW in our calcu- 
lations for small coupling strengths is questionable. 

Another way of understanding the formation of a CDW 
at quarter filling is to consider the energy gain by the 
Kondo effect. In Fig. HJ we show the kinetic energy, 
Ekin, and the energy of the exchange interaction, J(S-s), 
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Figure 3: (Color online) Upper panels: Transition tempera- 
ture over chemical potential for the CDW exactly at quarter 
filling. This plot illustrates the shrinking area of chemical 
potentials with increasing coupling strength, which stabilize 
the CDW at quarter filling. \i c is the chemical potential at 
the center of the CDW phase at quarter filling. Lower panel: 
Charge order parameter for different interaction strengths. 
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Figure 4: (Color online) Kinetic energy and J{S ■ s) (coupling 
energy between localized spin and conduction electrons) for 
two different interaction values J/W. By introducing a stag- 
gered potential for the sublattices A and B, we artificially cre- 
ate meta-stable solution with different occupation numbers, 
ua and ub- For all solutions the system is kept at quarter 
filling. 



for two coupling strengths calculated within DMFT at 
T = 0. Starting from the converged DMFT solutions 
for both interaction values, we artificially change the 
magnitude of the charge order by introducing a stag- 
gered chemical potential. This staggered chemical po- 
tential increases or decreases the difference in the oc- 
cupation numbers of sublattice A and B. For each ob- 
tained solution, we calculate the kinetic energy and the 
exchange energy. The calculated kinetic energies for both 



interaction values exhibit minima for ua = ns, corre- 
sponding to solutions without charge order. This can 
be expected, because neglecting the spin-electron cou- 
pling, the kinetic energy is supposed to drive the system 
into a paramagnetic state. On the other hand, the ex- 
change energy shows a very interesting behavior. For 
J/W = 0.4, which gives a CDW solution in Fig. [TJ 
J(S ■ s) decreases for increasing difference of the occu- 
pation numbers, ha — ns- Therefore, the exchange en- 
ergy favors for J/W = 0.4 a charge ordered state. The 
sum of kinetic energy and exchange energy results in a 
minimal energy for tia - hb ~ 0.2, which corresponds 
to the stable solution in our DMFT calculations with- 
out staggered potential. In contrast to J/W = 0.4, the 
exchange energy for J/W = 0.6 looks completely differ- 
ent. This energy has two local minima, one for a charge 
ordered state and the other for a state without charge 
order. The global minimum, and thus the energetically 
favored state, corresponds to the state without charge 
order. The sum of kinetic energy and exchange energy 
exhibits its global minimum for solutions without charge 
order, tia — tib = 0. 

Thus, wc arrive at the conclusion that a change in 
the behavior of the exchange energy, J(S ■ s), as a func- 
tion of the particle number is responsible for driving the 
system into a charge ordered state. We can qualita- 
tively explain this change by analyzing the Kondo en- 
ergy, E Kondo Wexp(-l/(Jp tF ))M The Kondo en- 
ergy depends via the DOS on the particle number. Thus, 
there might be parameter regions, where it is energeti- 
cally favorable to have different electron occupations on 
sublattice A and B in order to gain Kondo energy The 
importance of the Kondo effect and its dependence on 
the particle number is confirmed by the following two 
observations: DMFT calculations performed for the fer- 
romagnetically coupled Kondo lattice model, J < 0, at 
similar coupling strengths, do not show any sign of a 
CDW. Furthermore, DMFT calculations performed for 
the antifcrromagnctically coupled Kondo lattice, J > 0, 
but having a constant DOS, do also not show any sign 
of charge order. If the lattice has a constant DOS, the 
Kondo energy does not depend on the particle number, 
so that there is no energy gain by having different particle 
numbers on the two sublattices. 

In Fig. [5] we show spectral properties of the CDW 
at T = 0. The spectral functions of the CDW for both 
sublattices and for all coupling strengths J possess a gap 
at the Fermi energy, thus being insulating. This could be 
expected due to the long range order and commensurate 
filling. 

The qualitative shape of the spectrum can be under- 
stood in terms of a doped Kondo insulator. At half filling 
the Kondo effect splits the spectrum into two bands with 
a gap at the Fermi energy. For quarter filling, the Fermi 
energy is located at the center of the lower band. Because 
of the doubling of the unit cell due to the CDW, this 
lower band is again split into two bands with van Hove 
singularities at their band edges. This qualitative behav- 
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Figure 5: (Color online) Spectral function of the charge or- 
dered state for A- and B-sublattice for two different coupling 
strengths. The Fermi energy corresponds to cj = 0. 



ior can be explicitly seen for the J/W = 0.5 curve in 
Fig. [5] The upper band of the Kondo insulator, which is 
unaffected by the charge order, is located at w/W « 0.7. 
The gap of the Kondo insulating phase corresponds to 
the dip located at ui/W ~ 0.2. The lower band of the 
Kondo insulator, which is now split due to the CDW 
with van Hove singularities at the band edges, is located 
around the Fermi energy lo/W = 0. The charge order 
strongly affects the lower band of the Kondo insulator, 
which is filled differently on the two sublattices. For the 
smaller coupling strength J/W = 0.2, the gap of the 
Kondo insulator is strongly smeared out, which makes 
the identification of these structures more difficult. 



III. COMPETITION BETWEEN 
FERROMAGNETISM AND CHARGE ORDER 

As already noted in the beginning, most previous stud- 
ies found magnetism at quarter filling. Up to now we 
have completely suppressed magnetic ordering. In Fig. [5] 
we show the results of calculations simultaneously allow- 
ing for charge order and magnetism. In these calculations 
we allowed for any combination of magnetism and charge 
order on a two-site cluster. Thus, in principle also an- 
tifcrromagnctic solutions would have been possible, but 
only ferromagnetic results turned out to be stable in our 
calculations. The left panel and the middle panel show 
the strength of the spin polarization and the charge order 
parameter, respectively. From this we draw a qualitative 
phase diagram in the right panel. 

It should be noted that for these calculations we use 
mixing techniques within DMFT,— which help improv- 
ing the convergence towards the solution. Furthermore, 
there are parameter regions for which the DMFT calcu- 
lations without these mixing techniques at T = does 
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Figure 6: (Color online) Spin polarization (left panel) and charge order-parameter (middle panel) calculated from approximately 
300 data points for J/W = 0.35, simultaneously allowing for charge- and magnetic-order in two-site cluster. From these data the 
qualitative phase diagram (right panel) is drawn. The phase diagram comprises a pure ferromagnetic state at low temperatures, 
a transition region where charge and spin order can be found, a pure charge ordered state at intermediate temperatures, and 
a paramagnetic disordered state at high temperatures. 



not seem to converge, but get stuck in oscillating solu- 
tions. This could indicate that there are other competing 
solutions, which however could not be stabilized. 

From the panels showing the order parameters, it is 
clearly seen that the charge order parameter vanishes 
again at very low temperatures. The ground state is 
a magnetic (ferromagnetic) state without charge order. 
Thus, charge order and magnetic order are competing at 
T = 0. However, charge order can be found when increas- 
ing the temperature. The magnetic state vanishes at a 
lower temperature compared to the transition tempera- 
ture of the CDW. Therefore, exactly at quarter filling we 
will find a transition from a paramagnct to a charged or- 
dered state which is followed by a region of coexistence 
between charge order and magnetism, when decreasing 
the temperature. Eventually at even smaller tempera- 
tures, this coexistence region vanishes, and a magnetic 
state without charge order is formed. 

It should be noted that our calculations have been 
performed for an infinite dimensional Bethe lattice. Be- 
cause of the simplicity of the Bethe lattice, i.e. there 
are no closed loops, order patterns involving more than 
two sites and combining magnetic- and charge-order are 
unlikely. Furthermore, we have also checked different 
high-dimensional lattice geometries like the hyper-cubic- 
lattice, which exhibits a Gaussian DOS, and we found 
similar results as for the Bethe lattice. Thus, our results 
should be valid for bipartite lattices with similar DOS 
and large coordination number. We again want to stress 
that a lattice with constant DOS does not show charge 
order. Therefore, a curvature within the DOS is nec- 
essary. However, especially for low-dimensional lattices 
like the two dimensional square lattice, more complicated 
ground states allowing for a cooperation between charge 
order and magnetism might be favorable. In a recent 
study,— a combination of charge order and antiferromag- 
netism on four site clusters have been observed. Similar 
order patterns have also been observed for frustrated low 
dimensional lattices.— ~— However, the existence of these 
states strongly depend on the lattice geometry, and these 
states are not easily generalized to arbitrary lattices for 



higher dimensions, d > 2. 



IV. CONCLUSIONS 

We have shown that charge order can exist in the 
Kondo lattice model without any repulsive interaction 
between different lattice sites or between c- and /- elec- 
trons. For this purpose we have used DMFT, which is 
exact in infinite dimensions, but is supposed to be a good 
approximation for three dimensional systems. We have 
shown that a paramagnetic charge ordered state exists 
for small coupling strengths, J < 0.65 for a Bethe lat- 
tice, and vanishes as a first order phase transition for 
strong coupling. We have also shown that the charge or- 
dered state at quarter filling is an insulator, as can be 
expected due to the doubling of the unit cell. As an ex- 
planation for the existence of charge order, the non-linear 
dependence of the Kondo energy is proposed. 

Furthermore, we have studied the competi- 
tion/cooperation of charge order and magnetic order. In 
these calculations we have found that the ferromagnetic 
state without charge order is the stable state at zero 
temperature, and charge order can be only found for in- 
creasing temperature when the magnetic state vanishes. 
However, in all our calculations we have seen that the 
CDW has the higher transition temperature compared 
to the transition temperature of the magnetic state. 
Thus, the experimentally observed charge order can be 
explained without any repulsive inter-site interactions, 
as long as the system is in the Kondo regime (odd 
number of f-electrons), as it is assumed for a Kondo 
lattice model. 

Shortly before finishing this work, we became aware of 
a study by Misawa et al.— , who uses cluster methods 
and Variational Monte Carlo for analyzing charge order 
in the Kondo lattice model. 
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